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The effect of random surface roughness on quantum size effect in thin films is discussed. The 
conductivity of quantized metal films is analyzed for different types of experimentally identified 
correlation functions of surface inhomogeneities including the Gaussian, exponential, power-law 
correlators, and the correlators with a power law decay of the power density spectral function. The 
dependence of the conductivity a on the film thickness L, correlation radius of inhomogeneities 
R, and the fermion density is investigated. The goal is to help in extracting surface parameters 
from transport measurements and to determine the importance of the choice of the proper surface 
correlator for transport theory. A new type of size effect is predicted for quantized films with 
large correlation radius of random surface corrugation. The effect exists for inhomogeneities with 
Gaussian and exponential power spectrum; if the decay of power spectrum is slow, the films exhibit 
usual quantum size effect. The conductivity a exhibits well-pronounced oscillations as a function 
of the channel width L or the density of fermions, and large steps as a function of the correlation 
iy>2 ■ radius R. These oscillations and steps are explained and their positions identified. This phenomenon, 

fD ' which is reminiscent of magnetic breakthrough, can allow direct observation of the quantum size 

effect in conductivity of nano-scale metal films. The only region with a nearly universal behavior 
of transport is the region in which particle wavelength is close to the correlation radius of surface 
inhomogeneities . 

a' 

g 

o 

O I I. INTRODUCTION 

Progress in material technology, especially in nanofabrication, ultrathin film manufacturing, ultraclean and high 
^ • vacuum systems, etc., requires better understanding of boundary scattering in physical processes. The boundary 
effects should be an integral part of any study of quantum wires, wells, and films. Boundary scattering is especially 
important for transport in ultrathin and/or clean systems in which the particle mean free path is comparable to the 
' system size. 

Below we consider the effect of random surface roughness on quantum transport in quantized quasi-2£) systems 
' such as, for example, ultrathin metal films. The main issue is to find how sensitive is the transport along such film to 
the statistical properties of random surface inhomogeneities (thickness fluctuations) . An important by-product of our 
■ systematic comparison of different classes of random surface inhomogeneities is the prediction of a new type of size 
effect in quantized films. This effect manifests itself as large oscillations of conductivity a as a function of the film 
thickness L. In contrast to the usual quantum size effect (QSE), the peaks can be observed only at relatively large 
values of L. The distance between the peaks is large and is roughly proportional to L 2 . The observation of this QSE 
opens a new experimental method of identification of the type of surface roughness. 

The choice of quasi-2D systems is explained by a desire to avoid divergence of surface fluctuations and strong 
localization effects which are inherent to ID systems and make a systematic quantitative study of the effect of surface 
inhomogeneities on transport virtually impossible. In contrast to ID systems, the randomly fluctuating 2D surfaces 
are practically stable while the localization length in systems with weak surface roughness is exponentially large. 
(In general, the transport problems are more interesting in systems with weak rather than with strong roughness. 
Transport in systems with strong roughness is trivial: each wall collision completely dephases the particles and the 
mean free path cannot exceed the distance between the walls). 

The prevalent way to characterize the random surface roughness and/or thickness fluctuations is to use the corre- 
lation function of surface inhomogeneities 

C (s) = C (|s|) = (f(si)£(si + s)> = A- 1 1 £(si)£(bi + s)d Sl , (1) 

where the vector s gives the 2D coordinates along the surface, £(s) describes the deviation of the position of the 
surface in the point with 2D coordinates s from its average position, (£(s)) = 0, and A is the averaging area. Here 
it is assumed that the correlation properties of the surface do not depend on direction. Two main characteristics of 
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the surface correlation functions £ are the average amplitude (" height" ) £ and correlation radius (" size" ) R of surface 
inhomogeneities . 

Any transport theory for systems with rough boundaries should provide the explicit dependence of the particle 
mean free path (or the conductivity along the walls) on the correlator of surface inhomogeneities C (s). Without bulk 
scattering, the conductivity a is determined by the relation between three length scales: particle wavelength, A; the 
width of the channel, L; and the correlation radius of inhomogeneities, R. If the roughness is weak, the fourth length 
parameter, £, enters conductivity as a coefficient, 

°=—^f( A > L - R )- (2) 

Note, that this 2D conductivity differs by a length unit from the usual 3D conductivity and, as a result, has a 
dimensionality of conductance. 

The form of the surface correlator £ (s) can vary from surface to surface. Most of the theoretical calculations assume 
that this correlator is Gaussian. The numerical simulations, on the other hand, often rely on various generators for 
random rough surfaces without paying much attention to the correlation function of the generated inhomogeneities. 
Both approaches are not satisfactory since the experiments OBpSurface scattering and diffraction patterns show that 
real surfaces exhibit surface correlators of various formsua. Even one and the same film can exhibit various 
correlation properties on different stages of growth. As a result, the behavior of the functions / (A, L, R) in Eq. (||), 
which reflects the correlation properties of inhomogeneities, can vary from surface to surface even when the main 
correlation parameters £ and R remain the same. 

The correlation functions ([!]) are characterized by different long-range behavior that can be reliably identified in 
various surface diffraction and scattering experiments. What we would like to know is how sensitive is the particle 
transport to the form of the surface correlator. In contrast to surface diffraction and scattering data with angular 
and/or wavelength scanning, the transport coefficients arc integral parameters that include angular and wavelength 
averaging. This leaves the question of how sensitive is the conductivity to the shape of the surface correlator wide 
open. In addition, we are asking a question whether it is possible to identify the type of surface inhomogeneities 
from transport experiments in ultrathin films or multilayer systems without beforehand information on the form of 
the surface correlator. The interrelated question is, of course, to what extent one should pay attention to the details 
of the correlator of surface inhomogeneities in analytical or numericalr-transport calculations for particles with large 
mean free paths. The former issue has already been raised in RefsEHa for a small set of surface correlators on the 
basis of the Born approximation for wall scattering. Below we present a systematic study which is based on a more 
general transport formalism and involves a variety of classes of surface correlators. 

In short, we want to compare functions / (A, L, R) in Eq.(||) calculated for various types of the correlation functions 
£ (s) in a wide range of parameters. We start from degenerate ballistic fcrmions in quantized metal films. The choice 
is not arbitrary: transport in such systems involves the minimal degree of averaging (integration) and can be the most 
sensitive to the long-range properties of the surface correlators (Q) . 

Quantum size effect (QSE) in metal films is a, subject of intensive experimental study. Recent QSE experiments 
with quantized metal films include conductivity^, spectroscopy^, susceptibility!], and STM3 measurements. One of 
the signature features of QSE in metals is a pronounced saw-like dependence_jof conductivity on, for example, film 
thickness, <r(L). This dependence was predicted for both bulko and surface!!] scattering. Experimentally, QSE in 
conductivity was studied for metals in Refs.uiid (for earlier results see references therein). However, experiments on 
QSE in metals have to overcome a difficulty which one does not encounter in semiconductors. The period of the QSE 
oscillations in the dependence a (L) is usually small, almost atomic, 1/pf (below, except for final results, U = 1). For 
this reason typical experimental object are lead or semimetal films such as bismuth. Below we predict a new type of 
QSE with large-period oscillations of a (L) at relatively large values of pfL that could lead to observation of QSE in 
a wider group of metals. Large-period QSE oscillations have already been observed (see the second Ref.tl); however, 
sketchy experimental details do not allow one to identify reliably this observation as the new type of QSE predicted 
below. Our results can also resolve the long-standing controversy on the influence of the structure of the nanoscale 
film on its resistivity!! 

Recently, we developed a transparent semi-analytical formalism for transport in systems with rough boundaries that 
allows simple upifpipi calculations in a wide range of parameters and for variaus types of roughness with and without 
bulk scatteringEJii-l. This formalism unites approaches by Tesanovic et a/lla, Fishman and CaleckfiS, Kawabatal!3, 
Meyerov ich a nd S. Stepaniantsa, and Makarov et a£j (for a brief comparison between different theoretical approaches 
see Refs.t3E3). Below we apply this formalism with an explicit purpose of studying the dependence of the transport 
coefficients on the shape of the correlation function of surface inhomogeneities. The, well-defined limits of applicability 
of our approach to metal and semiconductor films are discussed in detail in RefsJ!3'0. 

Since the 2D mobility of particles is described by essentially the same equations as the exponent in the expression 
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for the localization length in films, our study provides the dependence of the localization length on the type of the 
correlation function of random surface inhomogeneities. 

The paper has the following structure. In the next Section we introduce various types of surface correlation 
functions. Section III briefly describes the transport equation used for conductivity (mobility) calculations in QSE 
conditions. The results of transport calculations for different types of correlators are given in Section IV. Conclusions 
and experimental implications are discussed in Section V. Appendix A contains useful analytical expressions for the 
power density spectral functions of inhomogeneities responsible for the behavior of scattering probabilities for different 
types of correlators. Appendix B deals with the positions of new type of QSE peaks. 

II. CORRELATION FUNCTION OF SURFACE INHOMOGENEITIES 

We consider an infinite 2D channel (or film) of the average thickness L with random rough boundaries 

x = L/2- 6(y, z), x = -L/2 + &(y, z). (3) 

(the walls are assumed hard with infinite potential). The inhomogeneities are small, £1,2 (y, z) <C L, and random with 
zero average, (£1) = (£2) = 0. Their correlation function Qk (s) and its Fourier image Cik (q), which is often called the 
power spectral density function or power spectrum, are defined as 

Cik (|s|) = (&(si)6c(si + s)) = A' 1 J 6(si)&(sx + s)dsi, (4) 



C«(|q|) = / d 2 s e^% lk (\s\) = 27r^ 



dk (s) Jo (qs) sds 



where s = (y, z) and q = (q y , q z ) are the 2D vectors. In homogeneous systems, the correlation function depends only on 
the distance between points |si — Sa| and not on coordinates themselves. The correlation functions Cu and C22 describe 
intrawall correlations of inhomogeneities, and £12 = (21 are the interwall correlations. Usually, but not always, the 
inhomogeneities on different walls are not correlated with each other, £12 = 0. Thus, everywhere, except for Section 
IVE, it is assumed that C12 = 0. To avoid parameter clutter, we also assume that the correlation parameters are the 
same on both walls, £11 = C22 = C Then the effective correlator contains 2£ (sLwith £ (s) given by equations below. 

Surface inhomogeneities exhibit a variety of types of the correlation functionsEm. To have a meaningful comparison, 
we consider the correlation functions that involve only two characteristic parameters, namely, the amplitude (average 
height) I and the correlation radius (average size) R of surface inhomogeneities. 

The most commonly used in theoretical applications is the Gaussian correlation function, 

C (s) = I 2 exp (~s 2 /2R 2 ) , ( (q) = 2tt£ 2 R 2 exp (~q 2 R 2 /2) (5) 

including its limit for small correlation radius R — ► 0, i.e., the 5-type correlations, 

C (s) = £ 2 R 2 5 (s) /s, C (?) = 2irfR 2 . (6) 

Sometimes, a better fit to experimental data on surface scattering is provided by the use of the exponential corre- 
lation function 

2tt/ 2 R 2 

C (.) = I 2 exp (s/R) ,C(Q)= — 3/2 , (7) 

(1 + q z R z ) ' 

or by the even more long-range, power-law correlators 



C (s) = 2 ^ rn , C (q) = 2^ 2 i? 2 (qR) (8) 

(i + s 2 /-R 2 ) 2/ j - 1 r(/x) MVy ' 



with different values of the parameter /z. The most commonly used are the Staras function with /u, = 1 and the 
correlator with fj, = 1/2 which has the exponential power spectrum £ (q), 

C (q) = 2it£ 2 R 2 exp (-qR) . (9) 

The use of the Lorentzian correlator, which differs from the definition (@) at /1 — > by the factor /1 in the numerator, 
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C (s) = 1 + S 2 /R 2 > C (<?) - ^ 2 R 2 Ko (qR) , (10) 

deserves a special comment. This correlator is often considered as " unphysical" . Its Fourier image ( |lC)| ) contains 
function Kq (qR) that diverges logarithmically at long wavelengths q — > 0. The issue to what extent the correlators 
are "physical" and can be reproduced experimentally is irrelevant in our context. For us, the fact that the Lorentzian 
correlator is sometimes used in calculations is sufficient enough to consider this correlator in the paper. To deal with 
the divergency, one can truncate the Lorentzian .correlator at large distances (the common practice is to make a cut-off 
at the distances about 0.1 of the system lengthfcl). Another option is to use the generalized power-law correlator (||) 
with small [i instead of the Lorentzian (|Io|). In order not to introduce additional parameters, we use the untruncated 
equation din) . Even though the divergence of Kq (qR — > 0) does not lead to any singularities in transport coefficients, 
the transport coefficients for Lorentzian surfaces (see below) often behave qualitatively different from systems with 
other types of random inhomogeneities, even from the systems (0) with small /i. (Sometimes, the divergence of the 
power spectrum £ (q) is associated with the fractal nature of the surfaceEl; to what extent our transport formalism can 
be used for films with fractal surfaces is an open question). 

The last class of correlation functions covers the power-law correlators in momentum space, 

C ( q ) = 2n£2R \^ , C (s) = i 2 } S/R)X K\ (s/R) . (11) 

(l + g 2 i? 2 ) 1+A 2 A L(1 + A) V 1 ' V ' 

ThepCorrelators from this group include the Lorentzian in momentum space A = that was observed in RcfS (see also 
Rcf.cl) and the exponential correlator (^ at A = 1/2. 

The constants in all these correlators are chosen in such a way that the value of £ (q = 0) = 2ir£ 2 R 2 is the same. 
This provides a reasonable basis of comparison for transport coefficients in films with all these different types of 
random surfaces. Indeed, the scattering cross-section for q — ► does not depend on the details of short- and mid- 
range structure of surface inhomogeneities. Therefore, at Fermi momenta qF — > (more precisely, at qpR <C 1), the 
transport coefficients should be the same for all random surfaces. (The only exception is the Lorentzian (|l0|) for which 
( (q) diverges at small q) . 

In what follows we compare the transport properties of the films - ([y]) in various ranges of the film thickness 
L, correlation radius R, and the particle wavelength Af — l/q F (or the 2D particle density N). 



III. TRANSPORT EQUATION FOR BALLISTIC DEGENERATE FERMIONS IN QUANTIZED FILMS 

QSE is caused by quantization of motion in the direction perpendicular to the film, p x — > 7r//L, and leads to a 
split of the energy spectrum e (p) into a set of minibands, e (p x , q) — * e (nj/L, q) = e, (q). For simplicity, we consider 
spherical Fermi surfaces ej (q) = ep, 



£3 (q) 



2m 



2me F - (irj/Ly 



1/2 



(12) 



. Qj = IF] 

where qj is the Fermi momentum for the miniband j. One can introduce the overall Fermi momentum as 

q F = 1/Af = (2me F ) 1/2 . (13) 

The relationship between this Fermi momentum q F and the 2D density of fermions N2 in quantized films is somewhat 
cumbersomeEj: 



N = J2 N 3 = (S/2tt) \q 2 F - (n/L) 2 (S + 1) (25 + 1) /6 
where S is the number of the occupied minibands, 

S = Int [gpL/Tr] 



(14) 



(15) 



If the density of fermions is the same as in the bulk, then N 2 = n^L where n 3 is the usual bulk density. Even in 
this case, the number of the occupied minibands S, according to Eqs. ( |l4| ) , (|l5|), is a complicated function of L. 
Asymptotically, at large S 



S = Int 



(3A 2 L 2 /vr) V3 



(16) 



4 



According to RefsJ3t2l, scattering by random surface inhomogeneities results in intra- and interband transitions 
e j (q) ~~ * e j' (l') with transition probabilities Wj.y (q, q') that are expressed explicitly via the surface correlation 
function £ (|q — q'|): 

(q,q') = Jj? [Cii + C 22 + 2C 12 f ?V f^V ■ (17) 
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The generalization to other, more complicated energy spectra is straightforwardlll. 
The transport equation for the distribution functions rij (q) , 



-i = ^ A Hj W 3 y[n f -n 3 ]5{e m -e^)^ (18) 
reduces, after standard transformations, to a set of linear equations 

3' 
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where rij — Vj5(e—eF)eE is the first angular harmonic of the distribution function rij (q) at q = qj, and 

Wjp 1 ^ {qj-Aj') are the zeroth and first harmonics of W {qj — q^) over the angle qTq^/- For some of the correla- 
tion function from the previous Section the angular harmonics can be calculated analytically (see Appendix A). For 
others, this calculation is performed numerically. 

The solution of Eqs. (|l9|) provides the conductivity of the film: 

e 2 x - 

3 

Equations (^) have simple analytical solution when the matrix rjy can be approximated by a diagonal matrix, 



T 33> 
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This happens when the matrix W^} is almost or exactly diagonal, w'J' ~ W^Sjji and 



2/mT j =*£W$-WW. (22) 

j' 



Then the conductivity rt21h is equal to 



e 2 ^ 2 2e 2 ^ 9 2 



,„ „, . TjQj = W^^y .,w {0 )-w iiy 

3 3 ^3' 33 3 

Such a diagonalization occurs in three physical situations. The simplest one is the one when only one miniband is 
occupied and 

— ^ ~ ~ 2 - ^ 1 (24) 



3h 2 m iyl Zh 2 m^ W [f-W^' 

The second case is the case of systems with large correlation length R^> L. In such systems the intraband scattering 
is much stronger than the interband one and the off-diagonal matrix elements Wjf are small in comparison with the 

diagonal ones (see Appendix A). Then both matrices wfy are almost diagonal, 
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W^^W^S jy , (25) 



and the expression for the conductivity ( |23| ) reads 

2e 2 \ - gj , , 

° ~ 3h 2 mi ^ W {0) - W (1) ' 

Such diagonalization of the matrices Wfy (^5|) at R~^> L can often be an oversimplification (see Section IV). 

The third situation with diagonal tJa, is the case of small qR. In this limit, the correlation function is a constant 
with the zero first harmonic, 



According to Eq.(|T^) . 



Wff = 2W {qR 0) , W\) ] = 



and 



a = 



2e 2 (L/irf ^ ( Lq } 



H 2S(S + 1) (2S + 1) c (0) ^ V 



Note, that all our surface correlators £ (s) are introduced in such a way that in the longwave limit £ (<7 - ► 0) they are, 
except for the Lorentzian (|l0|), equal to each other, £ (0) = 2ir£ 2 R 2 . This means that in this limit the conductivities 
(p8|) are the same irrespective of the shape of the correlator, 
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2e 2 1 (£ 2 /AR) 
ft 47r5(5 + l)(25+l)4-Uj/ ' 

(c/. Ref§). 

In all other situations Eqs.(|l9|) are not diagonal and should be solved numerically. 

The results for conductivity (mobility) also provide the exponent in the expression fpc the localization length 1Z 
that describes localization caused by particle scattering by random wall inhomogeneitiealj: 

n = Cexp[ranSD/h] (30) 

where C is the mean free path and the diffusion coefficient D is proportional to the conductivity a. 



IV. RESULTS AND DISCUSSION 



A. General comments 



As it is mentioned in Introduction, the 2D conductivity a of the film has the dimensionality of conductance and 
is described by a dimensionless function / in Eq. (^|) . This function, in turn, depends on the relation between three 
length scales - particle Fermi wavelength A^ = 1/<7f, the width of the channel L, and the correlation radius of the 
surface inhomogeneities R. The fourth length parameter, £, is perturbative and enters conductivity as a coefficient, 

<r = ^^f(A F ,L,R). (31) 

Note, that we consider only the contribution from surface roughness and disregard bulk scattering. As a result, the 
conductivity (^lj) diverges in the limit of vanishing inhomogeneities £ — > or R — > oo. The proper account of bulk 
scatterings eliminates this divergence. 

The dimensionless function f (Ap, L, R) depends only on the ratio of these three lengths. Of three ratios z = 
L/ Ap = q F L 1 x = R/ Ap = q F R, and y = R/L = x/z only two are independent, x = yz. Which two of these ratios 
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should be used as independent dimensionless variables depends on whether one wants to display the dependence of 
a on Ap, L, or R. The study of the dependence of the conductivity on film thickness, a (L), should be performed at 
constant and R. This means that a (L) is best displayed by the function fa (z, x) , 

2e 2 R 2 

a (L) = —-pft {z, x = const) , (32) 

for various values of x = R/Ap. 

Plots of the function fa (y) at constant values of z = qpL, 
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2e 2 L 

a(R) = —^f R (y,z = const) , (33) 

reflect the dependence a (R). Similarly, plots of the function fpf (z) at constant y = R/L, 

2e 2 L 2 

<r(qp) = -r--pfN{z,V = const) , (34) 

characterizes the dependence of conductivity on density of particles N or the Fermi momentum qp . 

Below we compare these dimensionless functions, //, (z) , fn (y) , and /jv (z) for various types of correlation functions 
in wide ranges of parameters. Needless to say, the results at x — > should coincide for all types of correlators except, 
maybe, for the Lorentzian, since, by design, all the correlation functions are the same in this limit [see Eq. (^)J. 

Curves in all figures below are labeled in a uniform way by the type of surface correlator used in calculations. 
Curves G correspond to Gaussian inhomogeneities (^|), Curves L describe the surfaces with Lorentzian correlations 
(|i~0|), Curves /jl±, /15, and u 9 give the results for the correlators (||) with u = 0.1; 0.5; 09, and Curves Ao, A5, and A 9 
correspond to Eq. ( |ll| ) with A = 0; 0.5; 0.9. Note, that correlator /15 has the exponential power spectrum (^) and that 
correlator A5 is actually the exponential correlator (Q). 



B. Dependence on the film thickness 

Figures 1 - 2 for the function f l (z, x = const), Eq. (|32"|), show the dependence of the conductivity a (V) for 
two different values of R/Ap, x = 1; 10, for various types of the correlation functions. The labeling of the curves 
G, L, Hi, /X5, //g, Ao, A5, Ag is explained in the end of previous subsection. The main feature of the curves, namely, 
their saw-like character, is well known. The sharp drops occur when the number of the occupied minibands, Eq. ( |15| ) 
changes by 1, i.e., in the points z = L/Ap = kir with integer k. The only unexpected feature is a "wrong" periodicity 
of the initial part of the Gaussian curve G at small values of z for x — 10 (see insert in Figure 2). This feature will 
be explained later. The Lorentzian curve L is different from others: at x = 10 the curve has already lost its QSE 
structure. 
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FIG. 1. Function / x (z,x — const), Eq.Jj^), at x — R/Af ~ 1 for various correlation functions. The labeling of the 
curves is explained in the end of Sec. IVA. Curve G: Gaussian correlator (g); curves /xi,/i5,/ig: power-law correlators (^) 
with = 0.1; 0.5, 0.9; curve L: Lorentzian correlator; curves Ao,A5,Ag: power-law correlators in momentum space ( pd| ) with 
A = 0; 0.5; 0.9 (A = 0.5 corresponds to the exponential correlator in the coordinate space (0)). The sharp drops occur when the 
number of the occupied minibands S, Eg. (Jig) , changes by 1, i.e., in the points z = L/Af — kn with integer k. 
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At these, relatively small values of x, the curves for all types of correlators have roughly the same shape though 
the exact values of the conductivity are different. (The curves /X5 and A5 are indistinguishable in both Figures 1 and 
2, and curves G and /Ltg are indistinguishable in Figure 1). To underscore this point, in Figures 3 and 4 we plotted 
instead of the curves fi, (z) the normalized curves /l (z) / fi, (z — z max ) with the normalization coefficients ensuring 
that the values of the normalized conductivity are equal to 1 at the highest values of z in the plot. Strikingly, for 
x = 1 (Figure 3) all the normalized curves with these 8 correlation functions lie within one bold line and are all 
indistinguishable with this resolution. For larger x, the difference is still insignificant: all the curves are compressed 
between curves G and L. The only anomaly is the loss of QSE structure by curve L. 




FIG. 2. The same as in Figure 1 for x — 10. The labeling of the curves is explained in the end of Sec. IVA. 
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FIG. 3. The same 8 functions /x (z, x = 1) as in Figure 1 normalized by their value at z = 110, }l (z) //l (110). All 
8 normalized curves are indistinguishable. The normalization coefficients are: curve G - Jl (HO) = 2.4 ■ 10 6 ; curve L - 
f L (110) = 1.39 ■ 10 6 ; curve /ui - ft (110) = 1.48 ■ 10 7 ; curve ^ - ft (110) = 3.61 ■ 10 6 ; curve fig - ft (110) = 2.42 ■ 10 B ; curve A 
- ft (110) = 2.69 ■ 10 B ; curve A 5 - ft (HO) = 3.65 ■ 10 6 ; curve A 9 - ft (110) = 4.54 ■ 10 6 . 

The main conclusion here is that the shape of the dependence a (L) at constant R and qF is not sensitive to and 
cannot provide any information on the type of the correlator at not very large values of R/ Ap. Since I is unknown 
and enters the conductivity as a coefficient, the absolute values of a (L) cannot serve as a clue either: experimental 
data on cr(L) at moderate R/Ap can be fitted by any type of the correlator by a choice of £. In this case, it is 
impossible to make any conclusion on the type of correlation function from transport measurements and it does not 
matter what correlator to use in theoretical calculations. Meaningful analysis requires some beforehand information 
on the correlation parameters. The only correlator that can be identified is the Lorentzian; however, this type of 
correlations is the least probable and might be " unphysical" . 

The situation changes dramatically at higher x = R/ Ap as is shown in Figures 5 [function (z,x = 400)] and 6 
[normalized function Jl (z, x — 400) //l (z = -z m ax, x = 400)] for the same 8 correlators (the labeling of the curves is 
explained in the end of Sec. IVA). 
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FIG. 4. The same 8 functions /l (z, x = 10) as in Figure 2 normalized by their value at z = 110, /l (z) / /l (110). All 8 
curves lie between normalized curves G and L and are barely distinguishable. The normalization coefficients are: curve G - 
f L (110) = 3.82 ■ 10 4 ; curve L - f L (110) = 1.17 • 10 4 ; curve fj,i - f L (110) = 1.48 ■ 10 5 ; curve ^ - f L (110) = 2.59 ■ 10 4 ; curve ^9 
-f L (110) = 1.32 ■ 10 4 ; curve A - f L (HO) = 6.95 • 10 3 ; curve A 5 - f L (110) = 2.61 ■ 10 4 ; curve A 9 - f L (HO) = 5.7 • 10 4 . 

We anticipated one feature, namely, the decrease in the amplitude of saw teeth with increasing x and even the 
disappearance of such teeth for the Gaussian correlator. The sharp drops in conductivity in the points where the 
number of the occupied minibands S increases by 1 is explained by opening of S new scattering channels associated 
with interband transitions in and from this newly opened miniband. Without the interband transitions, the increase 
of S by 1 results not in a sharp drop in a, but in an insignificant kink on the curve a (Li) as it is shown in the 
third of RefsJIa. The interband transitions are described by the off-diagonal components of the matrix of transition 
probabilities Wjf. With increasing R/Ap, these off-diagonal (interband) transition probabilities go to zero though 
with different rate for different types of the correlation function. The rate of decrease of the interband transition 
probabilities as a function of R/Af for different correlation functions is discussed in the Appendix. This rate is a 
good predictor for observing the saw-like shape of <j(L). The fastest decrease happens in the case of the Gaussian 
correlator; thus the curve for the Gaussian correlator should be the smoothest and should exhibit the smallest traces 
of the saw teeth. Therefore, the visibility of the saw teeth on the experimental curve can be a clue to the form of the 
correlation function. 
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FIG. 5. Functions /l (z,x = 400) for the same 8 types of correlators as in Figure 1. 

What is completely unexpected is the appearance of a new type of oscillation structure on a (L) in a limited range 
of z for the Gaussian and power-law correlators (curves G and fii in Figures 5,6). It looks as if there is a transition 
between two distinct regimes with several sharp oscillations in the transition range. The effect looks even more striking 
in Figure 6 for the normalized curves which, in contrast to Since Figure 5, is plotted in a linear scale. This new type 
of QSE requires an explanation. 

These new oscillations are not related to abrupt changes in the number of occupied minibands S (z): the oscillations 
are less sharp, have a much larger period, and, most important, appear only in a limited range of z where the number 







of occupied minibands S is already large. These new oscillations are observed for the correlators for which the 
interband transitions are the smallest and the saw-like structure is suppressed, namely for the Gaussian and power- 
law correlation functions. The power spectrum for these correlators £ (q) goes to zero exponentially at large q. Then 
one would expect that the off-diagonal (interband) transition probabilities are exponentially small in comparison with 
intraband scattering and that the conductivity can be well described by the " diagonal" approximation ( p6| ) that does 
not have an oscillation feature. This turns out not to be the case. 

The oscillations are indeed related to off-diagonal (interband) scattering probabilities Wjj>. A qualitative expla- 
nation of the effect and an estimate of the peak positions are the following. Scattering by surface inhomogeneities 
changes the tangential momentum by Aq ~ l/R. According to the momentum conservation law, this scatter- 
ing can cause the interband transition j j + 1 only when qj — qj+\ = Aq ~ l/R. If the miniband index 
j is relatively small and qj ~ 1/Ap, then qj — ~ (g| — q 2 + i) Ap/2. The energy conservation requires that 

q 2 — qj +1 = 7r 2 (j + l) 2 / L 2 — 7r 2 j 2 /L 2 ~ 2ir 2 j/L 2 . The combination of these conservation laws defines the peak 
positions Lj, which correspond to the opening of robust interband transitions j <-> j + 1 and which are given by 
equations L 2 ~ ir 2 jR\F. In dimensionless variables, this is equivalent to 

Zj ~ ITy^jx. (35) 

Accordingly, with increasing film thickness L the transition channel opens first for the electrons in the lowest miniband 
t\ (q) with j = 1. Note, that these are the grazing electrons which are responsible for the dominant contribution to 
the conductivity. Thus, the conductivity drops almost by half at the film thickness z-y ~ nx 1 ^ 2 where W\2 becomes 

1/2 

comparable to W\\ and the effective cross-section doubles. At higher value of L, z-i ~ 7r (2a;) , a new channel 
opens for the electrons from the next miniband j = 2 with p x = 2n / L and the conductivity drops again, and so on. 
The only difference is that the contribution of the electrons from the higher minibands falls rapidly with an increase 
in the band index j and the drops in conductivity <j{L), which are associated with the opening of new scattering 
channels for electrons from these minibands, become smaller and smaller. The number of the visible peaks on the 
curve a (L) and their relative heights give a good visual estimate of the number of "important" minibands and of 
their relative contribution to the conductivity. With further increase in the film thickness, when L becomes large, 
L > J!, the change of momentum Aq ~ l/R is sufficient to excite all interband transitions and the ordinary QSE 
with the saw teeth at the points z ~ nj is restored. 

The above explanation works for the films with the exponential decay of the power spectrum of inhomogeneities 
in which the size of inhomogeneities R is well-defined. In the films with a non-exponential power spectrum of 
inhomogeneities, i.e., with a more uniform distribution of inhomogeneities over the sizes R in momentum space, this 
new size effect cannot be observed because the particles from all minibands can always find the inhomogeneities of 
the right size that ensure the interband transitions irrespective of what is the separation between the walls. 
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FIG. 6. The same 8 functions /l (z,x = 400) as in Figure 5 normalized by their value at z = 110, /l (2) //l (110). The 
normalization coefficients are: curve G - /z, (110) = 1.84 • 10 5 ; curve L - /l (110) = 35.0; curve /ii - /l (110) = 8.78 • 10 4 ; curve 
/i 5 - h ( 110 ) = L2 5 ■ 10 4 ; curve /u 9 -/l (110) = 5.35 • 10 3 ; curve A - f L (110) = 3.16; curve A 5 - f L (HO) = 1.76 ■ 10 2 ; curve A 9 
- ft (HO) = 3.21 ■ 10 3 . 

More accurate explanation is the following. The off-diagonal elements Wjj> are functions of 
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Vjji = \qjR/h — qj/R/h\ = x 



^l~(irj/z) 2 -^l-(irf/z) 2 



and rapidly decrease with increasing Vjf, (see Appendix A). In general, the off-diagonal i/jj> is large at large R (or 
x) while the diagonal elements Vjj = 0. However, for large z (large S) some of the elements Vjji with small j, which 
are close to the main diagonal, could become small even for large x: 

2 

(j + 1 « z/n) ~?L£(2j + l) (36) 

(j changes from 1 to Int(z/7r)). Then at large z the transitions j «-> j + 1 can become noticeable and Eqs. ( |l9| ) 
become coupled. This coupling changes the solution of transport equation and, therefore, conductivity. According to 
Eqs. ( |T9| ) the coupling between the minibands j and j + 1 becomes noticeable, t^- +1 ~ rj- , when 

W® +1 (x, z) ~ wjf (x, z) - W$ (x, z) . (37) 

At fixed x, Eq. (|37]) can be considered as the equation for the values of z = Zj (x) at which one can observe the 
opening of transitions j <-> j ' + 1. The opening of such transition channels is accompanied by drops in conductivity. 
Since for the Gaussian and power-law correlators the interband transition probabilities Wjji depend exponentially on 
parameters Vjj> , these drops in conductivity are sharp and deep as illustrated in Figures 5,6. Solutions Zj (x) of Eqs. 
( |37| ) are discussed in Appendix B. At z = Z\ (a;), W\2 is the first of transition probabilities to acquire the "normal" 
order of magnitude. At z = z 2 (x), W23 becomes noticeable, then W34, etc. The amplitudes of the drops rapidly 
decrease with increasing j. In the end, when several interband channels with j <^ z/tt are open, a (L) becomes smooth, 
but with a much lower slope than in its initial part. The growth of transition probabilities for transitions j <-> j + 2 
does not result in new oscillations in a (L). In the points z (x) where Wjj+2 becomes large, W^ +2 ~ — Wjj , 
the states j and j + 2 are already strongly coupled via Wjj+i and Wj+i.j +2 . 

According to Appendix B, Eq.(^5[), the positions of the drops for films with Gaussian surface inhomogeneities are 
similar to Eq. (|35|): 

-1/4 



(i) w I V(2j + l)a; [in (xV2 (1 + 1/j))] • (38) 



The values zj (x = 400) = 33.4; 43.6; 51.8; 58.9; .... agree well with the positions of the conductivity drops on curve 1 
of Figures 5 and 6. 

For the surface with the power-law correlations of inhomogeneities @ the solution of Eq.(^7|) with logarithmic 
accuracy [Appendix B, Eq.(|60|)] again resembles Eq. (|35|): 



zj (x) = vrV(2j + l)x/4 i /, (39) 
v ~ In \x (1 + l/j) {2 In [a; (1 + 1/ j)]} m/2+1/4 ' . 

This expression is barely sensitive to \x. This almost complete independence of the peak positions from fi can be 
clearly seen in Figure 6. 

The difference between this new type of size effect and the usual saw-like QSE is dramatic. The saw-like drops in 
conductivity for usual QSE occur in the points z = kir with integer k and are direct consequence of quantization of 
momentum in thin films. The interband transitions are not germane to the existence or positions of this QSE and 
are responsible only for the amplitude of the conductivity oscillations. The drops in conductivity are equidistant with 
the period 7r along the z axis, i.e., are equidistant as a function of the film thickness. In contrast to this, the new 
QSE oscillations in Figures 5,6 are not related directly to the quantization of momentum and are a consequence of 
the exponential opening of interband transitions between minibands with small quantum numbers at certain values 
of the film thickness. The transitions in and out of higher minibands remain suppressed. (In some sense, the effect 
resembles magnetic breakthrough between separated parts of the Fermi surface in high magnetic fields). The peaks 
are roughly equidistant if plotted against z 2 ; weak deviation from periodicity is due to a logarithmic terms in Eqs. 
( p8| ) , ([39]). The period of the new QSE is much larger than for the usual QSE. The large period of oscillations can 
open the way to direct observation of QSE in transport measurements in metal films in which usual QSE has atomic 
period and can hardly be observed. There is a strong possibility that the conductivity oscillations reported in the last 
of Refs. ■ are actually this new type of QSE. 

The initial part of the curves G, fii in Figures 5,6 for a (L) is described analytically by Eq.(j2(|) with appropriate 
values of W from Appendix A. This curve is close to the power law a oc i( 5 + Q ) (small a depends on x) and to 
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experimental data of the third Ref. ? . After the region of new QSE oscillations, the curves are again smooth, but with 
a much smaller tangent. We do not have an analytical description for this regime. The numerical approximation 
can be done equally well by either a = A + B ■ L 1+l3 with small also depends on x) or a quadratic expression 
a + b ■ L + c ■ L 2 . This behavior explains the experimental datao and the last Ref.Q. As a result, the power-law 
dependence of a (L) is qualitatively different for ultrathin and imp; -thicker films. This type of behavior is different 
from the earlier studied behavior of a (L) at small x = qpR l.Eni3'E3. 
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FIG. 7. Functions $l (z, x — const) for Gaussian correlation of surface inhomogeneities normalized by their value at z = 157, 
/l {z) J fh (157). The values of x and normalization coefficients are: curve 1 - x = 1, / (157) = 6.9 ■ 10 6 ; curve 2 - x = 10, 
/ (157) = 9.9-10 4 ;curve3-a: = 25, / (157) = 4.6 ■ 10 4 ;curve 4 - x = 55, / (157) = 3.8- 10 4 ; curve 5 - x = 100, /(157) = 4.75 TO 4 ; 
curve 6 - x = 200, / (157) = 9.1 ■ 10 4 ; curve 7 - x = 400, / (157) = 2.3 ■ 10 5 . 

The initiation of this new type of oscillations with a large period can be seen on the initial part of curve G in 
Figure 2 for x = 10. With growing z these new oscillations get overtaken by the standard QSE. The transition from 
standard to new QSE is illustrated in Figure 7 that contains normalized "curves G" for the Gaussian inhomogeneities, 
j 'l(z,x — const) I Jl{z — 157), for x = 1; 10; 25; 55; 100; 200; 400. It is clear from these curves how usual QSE is 
replaced by new oscillations with increasing x. The "transitional" curve for x — 55 is especially interesting: it shows 
new QSE at smaller z and a restoration of standard QSE at higher z. This restoration occurs when a noticeable 
number of interband transitions become open at higher z. It seems that such restoration does not happen on curves 
x > 50. This impression is wrong. Such restoration indeed occurs for curves x = 100; 200; 400, but at values of z 
that are much larger than those in the Figure. At very large x, all curves /l (z,x = const) consist of four parts: 
rapid increase at small z, region of new QSE oscillations, smooth monotonic part, and the region of relatively smooth 
standard QSE oscillations at the largest values of z. With increasing x, the amplitude of new QSE oscillations and 
the length of region separating new and old QSE increase rapidly. 

C. Dependence on the correlation radius 

The dependence of the conductivity on the correlation radius of surface inhomogeneities, a (_R), is best illustrated 
by the function fa (y, z = const), Eq. (|35|). Since the number of the occupied minibands S does not depend on the 
correlation radius of inhomogeneities, the curves fa (y) at constant z do not exhibit the saw-like structure. Instead, 
the two main features are the presence of the minimum in fn (y) and the step-like structure that corresponds to the 
oscillations in Figures 5,6. 
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FIG. 8. Function fa (y,z = 64.4), Eq.(|33|), near the minimum at yz ~ 1 for various surface correlators. The labeling of the 
curves is explained in the end of Sec. IVA. 

The scattering of fermions by surface inhomogeneities is most effective at R/ 1 Ap ~ 1, i.e., at y ~ 1/z. This leads to a 
minimum of the conductivity a (R) at such values of y. At R/ Ap <C 1 the particle wavelength is much larger than the 
size of surface inhomogeneities and the scattering is almost specular and does not contribute to the formation of the 
mean free path. In the opposite limit R/Ap 3> 1 the walls are flat on the particle length scale and surface scattering 
also does not limit the effective mean free path. Therefore, at z = const the conductivity a (R) for non-divergent 
correlators is infinite in both limits y — > and j-too with a minimum around y ~ 1/z. The curves fa (y) close to 
this minimum are plotted for different correlators in Figure 8 (z = 64.4; the labeling of the curves is explained in the 
end of Sec. IVA). It is important that the position of the minimum, its width, and even the order of magnitude of 
the function fji (y) in the minimum are roughly the same for all types of surface correlators. This is, probably, the 



most universal feature of the system. The only correlator that does not display a well-defined minimum is (11) with 
A = (the Lorentzian in momentum space; curve Ao). This feature is related to the logarithmic divergence of this 
correlator in "real" space. This feature is especially interesting because the surfaces with such inhomogeneities were 
observed in experimento. 

The drops in a (L) at large z — Zj (x), which are analyzed in the previous Section (Figures 5,6), correspond to the 
points yj (z) on the curves fu (y). The positions of these points yj (z) are implicitly determined by Eq. ( |3q ) and ( |39| ) 
for the Gaussian and power-law correlations provided that x = yz. These values of y are far away to the right from 
the minimum in the curves a (R) and cannot be presented in the same figures. The feature that corresponds to the 
oscillations from the previous Section is clearly seen as a set of steps in Figure 9 for the same value of z as in Figure 
8, z = 64.4 on curves G and ^5 for Gaussian and power-law inhomogeneities. For the surfaces with the Gaussian 
inhomogeneities , the first interband transition W\2 becomes visible for z — 64.4 at j/i ~ 25, the next one at yi ~ 14, 
and so on. At these values of y one can see well-pronounced steps on the curve G in Figure 9. The same feature, 
though barely discernible, is also observed for the power- law correlator /15. 

For comparison, curves L, Ao, and A5 do not exhibit any anomalies. Interestingly, the curve for the Lorentzian 
inhomogeneities is the only one that decreases with increasing y after the initial increase at small y (Figure 8). How 
is this feature related to the peculiarities of the Lorentzian that have been discussed in Section II is unclear. The 
curve Aq remains essentially flat. 
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D. Dependence on the Fermi momentum and density of fermions 



The dependence of the conductivity a on the density of fermions N or their Fermi momentum qp is best displayed 
by the function f N (z) at constant y = R/L, see Eq. |||). This dependence a (N) is similar to a (L). The function 
a (N) exhibits a clear saw-like structure of usual QSE at not very high y for all correlators. With increasing y, the 
saw teeth disappear first for the Gaussian correlator G, then for the power-law correlators but persist for the 
power-law correlators in momentum space A^. Instead, at large y the functions /jy (z,y = const) for Gaussian and 
power-law inhomogeneities exhibit new type of QSE oscillations similar to that for Jl (z, x — const) in Sec. IVB. The 
positions of these oscillations can be found from Eqs.(|38|) , ( j39[ ) after substitution x = yz. 
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FIG. 10. Normalized function /at (z,y = 20) Eq.(^iJ), /at (z) / /at (z = 126), for three surface correlators. The normalization 
coefficients are: curve G - f L (126) = 1.1 ■ 10 9 ; curve ^ 5 - fh (126) = 4.5 ■ 10 7 ; curve A 5 - f L (126) = 1.4 ■ 10 4 . 

This effect is illustrated in Figure 10 (the labeling of the curves is explained in the end of Sec. IVA). The figure 
presents functions fx (z, y = 20), Eq. (]34|) , for the Gaussian (curve G) and power-law (fi — 0.5; curve /i 5 ) correlators, 
and for the correlator with a power-law power spectrum (A = 0.5; curve A5). To compensate for different orders of mag- 
nitude of the data for these correlators, the functions are normalized by their values at z = 126, Jn (z) / Jn {z — 126). 
Curve A5 exhibits a saw-like behavior typical to the usual QSE with period 7r along the z-axis. Curves G and ^(5 
exhibit new QSE oscillations with a much larger period. 
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E. Interwall correlation of inhomogeneities and quantum size effect 



Surprisingly, the possibility of interwall correlation of surface inhomogeneities gives an interesting insight into usual 
and new QSE and provides an additional proof for our explanation of QSE oscillations reported above. The study of 
the effect of interwall correlation of inhomogeneities has been initiated in Ref.E3 for Gaussian correlations. Below we 
supplement those results for other types of surface correlators with-an emphasis on new QSE. 

To decrease the number of parameters, we assume that, as in Ref.E3, the correlation functions of inhomogeneities on 
both walls Cii and (22 are given by the same function, (n (s) = £22 (s) = C ( s )- The structure of the interwall correlator 
of inhomogeneities £12 (s) is assumed to be the same as for the intrawall correlations with the same correlation radius 
R. However, the amplitude a of the interwall correlations is different from the intrawall ones, 

Cu = C22 = C(a), Cm(«) =<(s). (40) 

To compare the effect of such interwall correlations for different classes of the function £ (s), we calculate the relative 
change of conductivity a (i.e., functions /l, fa, /at) caused by introduction of such correlations 

t ^ f (a) - f 

«/>(*) = (41) 

where and / are the functions Jl,r,n calculated with and without interwall correlations. An additional benefit 
is that the functions for all types of correlators are automatically normalized thus eliminating a difference by 
orders of magnitude between the functions /l,r,n for different types of correlation functions. 




FIG. 11. Relative change <f>f n) , Eq .(fill) , of the function /l (z,x — 1) , Eq. (^i|) , for the interwall correlation amplitude 
( fIo| ) a = 0.75 for various correlation functions of surface inhomogeneities. The labeling of the curves is explained in the end of 
Sec. IVA. All the curves exhibit almost identical oscillations as it should be for well-developed usual QSE. 

In the presence.of such interwall correlations, the transition probabilities Wjf (q, q') ( p7| ) become proportional, in 
accordance witblij, to 

l + a(-l) J+/ ]C(|q,-4|). (42) 

The most interesting effects of the interwall correlations are related to the oscillating structure of the term with a in 
Eq. (H|). If the interband transition probabilities W^y (q, q') are large, i.e, if ( (|qj — q^-, |) is not small for j' ^ j, 
then the contribution of the term with a in ( fl2"| ) has a different sign for different Wj? depending on whether j + j' 
is even or odd. This should result in an oscillating structure of the function <^>W (flj) as a function of the number of 
occunied nhnibands 5, i.e., as a function of film thickness L (the existence of such oscillations was first reported in 
Ref.E3 for Gaussian inhomogeneities). The period of such oscillations should be equal to that for standard QSE and 
their amplitude should decrease rapidly with increasing L. Since our explanation of the standard QSE ties it to large 
interband transitions, the oscillation nature of the function <j)^ (|dj) should exist in the same range of parameters as the 
standard QSE. In accordance with Sec. IVB, these oscillations should be noticeable for the function <fi^ (z, x — const) 
at small x for all types of surface correlators. This is illustrated in Figure 11 (x — 1) for the correlators G, L, A5, 
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/i5- The figure is plotted for a = 0.75. The similarity of the functions <f> L ' (z, x = 1) is striking, but not surprising. 
The flat part of all curves at small z is explained below. At higher values of x, the interband transitions (off-diagonal 
Wjji become more and more suppressed. When the interband transitions become negligible, the only non-zero 

scattering probabilities are diagonal Wjj that are proportional to 2 [1 + a] C (|qj — q^j) Eq.([42|). Since all Wjj are 

scaled by the same factor 1 + a and the conductivity is inversely proportional to W , the function <j>^ (z) in the absence 
of the interband transition becomes a constant, 

If a = 0.75, the value of this constant is 0^ ' 75 ' ) (z) — —3/7. Eq.([43"|) also describes the initial part of all curves <pf/ ( z ) 
for all values of x at small z when only the first miniband is occupied, S = j = j' = 1. This explains all curves in 
Figure 11 have identical flat pat at small z. 
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FIG. 12. Relative change (j> ( °' 75) , Eq.(|lJ) , of the function f L (z,x = 400) (§|) for the interwall correlation amplitude @ 
a = 0.75 for four correlation functions of surface inhomogeneities. The labeling of the curves is explained in the end of Sec. 
IVA. Curve As exhibits oscillations in accordance with usual QSE for curve As in Figure 6. Curves G, L, and fis are flat at 
small z, 4> { l 75) = -3/7, Eq.©. Oscillations on curves G and ^ confirm the explanation of the new QSE as an exponential 
appearance of transitions j < > j + 1 at certain values of z. 

Figure 12 illustrates <jyf^ (z,x — const) at x = 400 and a = 0.75 for several correlators. At this value of x, the 
exponential correlator A5 (pi) exhibits, according to the results and explanation of Sec. IVB, the usual QSE. Therefore, 

the function </>^°' 75 ' > (z, x = 400) for this correlator should have an oscillation structure; this is clearly seen in Figure 12. 
The Gaussian and power-law correlators G and a 5 , according to Sec. IVB, ensure the absence of interband transitions 
at small and moderate z where the function (fi^' 75 ^ = —3/7 in Figure 12. Our explanation for the new type of QSE in 
Sec. IVB is an abrupt sequential appearance of noticeable interband transitions W12, W23, W34,, etc. at certain values 
of z = Zj. Since for the term with a in Eq.(^2|) is negative for all transitions f = j ± 1, one should observe spikes in 

conductivity and, therefore, in the function 4>l \ at 2 = Zj. In some sense, Figure 12 provides the best illustration for 
our explanation of new QSE. 

Figure 12 also provides an insight into anomalous behavior of conductivity for Lorentzian correlation of inhomo- 
geneities (|hJ), curve L. At z < 30, the interband transitions are suppressed and (fi'f' 75 ^ = —3/7. At higher z, the 
interband transitions become more noticeable and start increasing, but very slowly. Why does the curve remain 
smooth when a sufficient number of transitions is already visible, is still a puzzle. A possible explanation is that 
oscillations should appear only at very large S (or z) when their amplitude should be vanishingly small. 

V. SUMMARY AND CONCLUSIONS 

In summary, we compared the behavior of conductivity for various types of surface correlators in a wide range 
of parameters. The following conclusions can be important when analyzing the experimental data or discussing 
theoretical predictions. 
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The rough shapes of the curves of the transport coefficients are similar at small and moderate R for all types of 
correlators though the orders of magnitude of the transport coefficients and more fine details of the curves can 
be different. To make any definite conclusions from the rough shapes of the experimental curves, one should 
have at least some idea of the type of the correlation function of surface inhomogcncitics and/or the value of 
the correlation radius R and the amplitude of inhomogeneities I. Since t plays the role of a scaling parameter, 
getting the values of parameters of surface inhomogeneities from experimental data on transport without any 
additional information on the correlation of inhomogeneities could result in mistakes by orders of magnitude. 
In the same way, the use of the wrong correlator in theoretical calculations could result in absolutely wrong 
predictions without evoking any warning signals from comparison of the rough shapes of experimental and 
theoretical curves. 

The most universal feature is the shape of the curves and order of magnitude of a (R) near the minimum at 
R/Af ~ 1. This minimum allows experimental evaluation of the correlation length of surface inhomogeneities 
R without any assumptions about the type of the correlation function. 

The shape of the curves &(L), u(N), and a (R) becomes very sensitive to the type of surface correlator at 
large correlationjadius of inhomogeneities R. Experimentally, this is important for better quality films (see, for 
example, in Ref.E3) in which STM and other usual methods are not well-suited for the study of the long-range 
behavior of the thickness fluctuations. Here transport measurements can be used as a good alternative for 
identification and analysis of the thickness fluctuations. 

The underlying reason is very high sensitivity of coupling between quantum well states with low quantum 
numbers to film thickness and the long-range behavior of the thickness fluctuations. This phenomenon is quite 
general and should lead to observable effects not-only in metal films, but for other types of quantum wells such 
as semiconductor films or quantum wave guidesu. 

The persistence of the saw-like dependence of the transport coefficients on the thickness of the film, Fermi mo- 
mentum, or the density of fermions should signal the long range nature of the surface correlations in momentum 
space C (?)• The observation of the saw- like structure for R > L is a distinct signature of the power-law decay 
of the power spectral density function ( (q) though, by itself, is insufficient to make conclusions about the index 
in this power law. The easy suppression of the saw-like behavior points at the exponential decay of the power 
spectral density. The rate of this suppression is significantly different for simple exponential and Gaussian decays 
ofC(ff). 

Thickness fluctuations with Gaussian correlations and correlations with exponential power spectrum lead to a 
new type of QSE in a (L), a (N), and a (R) for surface inhomogeneities of a relatively large size R. This new 
QSE produces large oscillations in a (L) and a (N) and steps in dependence a (R). The spacing between these 
new QSE anomalies provides important direct information on the correlation parameters of inhomogeneities. 
The peaks are almost equidistant if plotted against z 2 . 

In contrast to the usual saw-like QSE, the new QSE oscillations are not related directly to the quantization 
of momentum and are a consequence of the exponential opening of interband transitions between minibands 
with small quantum numbers at certain values of the film thickness. In some sense, the effect is reminiscent 
of magnetic breakthrough that describes the opening of transitions between disconnected parts of the Fermi 
surface. 

Large period of new QSE oscillations opens the way to direct observation of QSE in conductivity of quantized 
metal films and may be responsible for experimental data in the second Ref.13. An additional experimental 
signature should be the appearance of these new QSE oscillations only at relatively large values of the thickness 
of quantized metal films. 

The Gaussian correlation of inhomogeneities affects particle transport in a unique way. First, the values of 
the transport coefficient are, except for the smallest correlation radii, larger than for other, slower correlators 
by orders of magnitude. This is explained by this correlator having the shortest tails resulting in the least 
effective scattering. Second, this type of correlation does not exhibit a saw-like dependence of the transport 
coefficients on the system parameters except for small correlation radii R. Third, this type of correlation of 
the surface inhomogeneities leads to the above-mentioned new type of large-scale oscillations of the transport 
coefficients. The combination of these features can make the Gaussian correlator readily identifiable in transport 
experiments. 
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• The Lorentzian correlation of inhomogeneities in configuration space is also readily identifiable by several ab- 
normal features. The combination of these features could be another manifestation of an "unphysical" nature 
of this correlator. If possible, this correlator should be avoided in theoretical and computational models. A 
power-law correlator (|) with small index [i can serve as a good replacement in the calculations. 

• The results explain the observed difference in power-law regimes of the thickness dependence of the conductivity 
a (L) between ultrathin and more thicker films. 

• The relative contribution of the interwall correlation of surface inhomogeneities strongly depends on the type of 
QSE. For usual QSE, the contribution of the interwall correlations is a rapidly decaying oscillation function of 
the film thickness. For QSE of the new type, this contribution is constant in a wide range of small and moderate 
thicknesses, and becomes an oscillating function with a big period in a limited range of large thicknesses. 
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VII. APPENDIX A. TRANSITION PROBABILITIES 

Various correlation functions from Section II allow different degrees of analytical calculations of the scattering 
probabilities. The angular harmonics of the correlation function £ (|q — q'|) in the transport equation ( JTgj ) are defined 

as 



1 00 
C (|q - q'l) = ^C (0) («, q') + E ^ («« <*') cos (s X ) ■ 



(44) 



c 



C(x) cos (sx) dx 



where x is the angle between the 2D vectors q and q'. 
The harmonics for the Gaussian correlator (^|) are 



C (1) (<&,«&') = ^l 2 R 2 



-W/o (QQ') 
-W'h {QQ 1 ) 



e -(Q-Q')7 2 , 
£ -(Q-Q') 2 /2 



(45) 



where Q — QjR, Q' — qj'R- Note, that in RefsJla we used equivalent expressions with hypergeometric functions 
instead of modified Bessel functions. Expressions in square brackets in Eqs. ( f44| ) are smooth functions of Q and 

Q'. The exponential coefficients, exp — (Q — Q') 2 /2 , on the other hand, are rapidly going to zero for large qR if 

1j ^ Qj'- This explains why the off-diagonal scattering probabilities Wjj> are much smaller than the diagonal ones 
at large qR. Such a drastic difference between interband and intraband scattering probabilities is a unique feature of 
the Gaussian correlator. The physical consequences are discussed in Section IV. 
For the exponential correlator (Q) the harmonics are 



" 2 R 2 E(n) 



1 + iQ-Q'Y ^1 + {Q + Q') 



(46) 



t 2 R 2 I 1 



Q' 2 )E(n)- (l + (Q-Q'y)K(Q,) 



QQ' 



l + (Q-Q0 \n + (Q + Q') 



C (0) (qj,Qf) 

C (1) (qj,q?) 

n = 2jQQ'/ [1 + (Q + Q* 

where E and K are complete elliptic integrals. Here the diagonal and off-diagonal transition probabilities (probabilities 
of the intraband and interband scattering) differ mainly by the terms l+(Q — Q') in denominator that are insignificant 
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in comparison with the exponential factors for the Gaussian correlator above. The physical consequences are discussed 
in Section IV. 

The power-law m) correlation functions correspond to 



'(0) 



Kfj,-\-m (Qmax) ^fi-\-m (Qmin) 



jmax 



Q 



min 



(47) 



/ (cos <f>) [Q 2 + Q' 2 - 2QQ' cos <f\ M d<t> 
Jo 

K fi+m (Qmax) ^Lt+m (Qmin) 



C (i) = (/i + m) 



m=0 



; max 



min 



/•27T 

/ C£ (cos <t>) [Q 2 + Q' 2 - 2QQ' cos 0] M cos < 
Jo 



where are the ultraspherical (Gegenbauer) polynomials, and Q max = max(Q, Q') and Q m ; n = 

mm(Q,Q'). The off-diagonal transition probabilities disappear exponentially at large \Q — Q'\, approximately as 
(\Q — Q'ly 1-1 ^ 2 ex P ( — \Q — Q'\): i- e -; much slower than for the Gaussian correlator ( [45|) but faster than for the cor- 
relator @. 

The integrals in Eqs.(0) can be simplified for the Lorentzian correlator: 



C (0) (Q, Q') = &Trt 2 RKa (Q max ) I (Q min ) 
C«(Q, Q') = ^l 2 RK x (Q max ) h (Qmin) 



(48) 



Note, that the function Kq (Q) diverges logarithmically at Q — * 0. This divergence is discussed in Sections II and IV. 

The expressions for the harmonics ( f47j ) can also be simplified for the Staras correlator, fj, = 1, when C* (cos</>) = 
sin [{n + 1) 4>] I sin^, 



I C} n (cos 4>) d<p = [0, m = 2k + 1; 2vr, to = 2fc] , 
Jo 

p2rr 

/ (cos </>) cos ^ d</f> = [0, to = 2fc; 27r, to = 2fc + 1] , 
Jo 



and the harmonics (|47|) reduce to the rapidly converging sums of the Bessel functions with alternating coefficients. 
For all other power-law correlators with different values of fi the integration should be performed numerically. 

The last group of correlators involves power-law behavior in momentum space, Eq. (|ll|) . This group includes the 
Lorentzian in momentum space A = that was observed in Ref.0 and the exponential correlator (|7|) , I 
In general, the angular harmonics are 



£(0) 



4tt£ 2 R 2 



1 + {Q 2 -Q' 2 ) 2 + 2{Q 2 + Q' 2 ) 

Anl 2 R 2 /\ 

1 + (Q 2 - Q' 2 f + 2 (Q 2 + Q' 2 ) 



(l+A)/2 



Pa (ft) 



at A = 1/2. 
(49) 



(l+A)/2 



Pa 



n = (1 + Q 2 + Q' 2 ) /yjl + (Q 2 - Q' 2 ) 2 + 2 (Q 2 + Q' 2 ) 

where P" (fi) are the associated Legendre functions of the first kind. Note, that the argument fl of the Legendre 
functions in our expressions is larger than 1. One should be cautious when doing calculations with expressions (^9[): 
some of the handbooks (and software packages, e.g., Mathematica) do not use the same normalization for Legendre 
polynomials and Legendre functions, i.e., for functions P" (f2) with integer and non-integer A. 
In the case of the Lorentzian in momentum space, A = 0, 



Ciit (*) = e 2 K (s/r) , c(Q) 



2^ 2 P 2 
l + (Qi?) 2 ' 



(50) 



the harmonics 
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4tt£ 2 R 2 



C (0) {w) 
C (1) fo.fcO 



i + Q 2 + Q 12 + Vi + (Q 2 - Q' 2 ) + 2 {Q 2 + Q' 2 ^ 

Note, that this correlator diverges in real space at s — > 0. 



1 + (Q 2 -Q' 2 ) 2 + 2(Q 2 + Q' 2 ) 
8ir£ 2 R 2 QQ' 

1 + (Q 2 - Q' 2 ) 2 + 2 (Q 2 + Q' 2 ) 
1 



(51) 



VIII. APPENDIX B. POSITIONS OF NEW QSE OSCILLATIONS. 

The peak positions are determined by the condition that the absolute value of the diagonal and the first off-diagonal 
matrix elements in transport equation (^) become comparable: 



l/r : 



1M. 



Rewriting this condition via transition probabilities Wjp 1 ^ (q, q') we get 

W$ (x, z) - W$ (x, z)] + W$ (x, z) ~ W$ +1 (x, z) 

]'¥=] 



(52) 



where W^i (qj,qj') are the zeroth and first harmonics of W (q^— q^-,) over the angle qjq^ that can be expressed 
explicitly via the surface correlation functions [see Eq. ([[t]) and Appendix A] . For large qjR, the off-diagonal scattering 
probabilities Wjj> are exponentially suppressed for Gaussian and power-law inhomogeneities, Eqs. J4E| ) ,(p7|) : ~ 
> WfL ~ WfV+i. With a logarithmic accuracy, the condition ( |52"| ) corresponds to the equation 



Taking into consideration the asymptotic behavior for modified Bessel function in Eq. 
lator, Eq.(p3]) can be reduced to 



f _ (i + if 

2Qj y/QjQ 



exp 



1 s 



(53) 

for the Gaussian corre- 
(54) 



where Qj = x\J 1 — (irj/z) . When z/nj 3> 1, we can put Qj « Qj+\ ~ x in the denominator. The exponent should 
be evaluated more carefully: Qj — Qj+i ~ xir 2 (2j + 1) /2z 2 . Then Eq. (El) yields the following values of the peak 
positions: 



z 3 ( x ) = 77 



In 



Vs(2j + 1) 

1/4- 



(55) 



Since x = yz, these peak positions Zj (x) can be also used to get the peak positions for the conductivity at fixed y, 
Zj (y) as a solution of the following algebraic equation: 



Zj (y) 



y(2i + i) 



In 



1/2- 



(56) 



Similar but more cumbersome calculations, can be performed for the power-law correlators (||) . For example, if 
M = 1/2, Eq. d||) reads 
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ir/2 tt/2 

4? 2 J exp (-2Qj sin*) sin 2 t eft = 2 (j + l) 2 ^ exp (^-\Jv 2 + ^QjQj+i sin 2 eft, 



(57) 



where we introduced i/ = Vjj+i = Qj — Qj+i. For large Qj, an asymptotic estimate for the integral in the left-side is 
Rough asymptotic estimate for the integral in the right-side of the equation is 



exp -\v 2 +4Q,Q,+it 2 



dt 



VT^t 2 2^Q 



2 \/ QjQj+i 

===== / exp [-\/v 2 + y 2 ) dy 
jQj+i J J 



In order to estimate this integral, we can substitute \/v 2 + y 2 by 



V" 2 + : 



V, for y < v 
y, for y > v. 



Then 



OO 

4 / cxp i-^^) d y * <?+ !) • 



This leads to the following estimate for peak positions 



Zj (X) = 7T 1 



lx(2j + l) 



2vi 



where Uj is the root of the transcendental equation 

Vj = 2 In Aj + In (1 + Uj) , , Aj = x (1+ 1/j) . 
The last equation can be solved by iterations: 



2 In A, = In [2 In A + 1] 



Finally with a logarithmic accuracy, the solution of Eq. ([57]) for the positions of peaks becomes 



x(2j + l) 



1 \ In [xVln (x(l + (1 + 1/j) 
Similar asymptotic estimates for the power-law correlators with arbitrary \i yield 



ar(2j + l) 



\ In |Aj (2 In ^ 



W2+1/4 



3 j yr^ + 5/2) 

It is clear from Eq. (|6^) that the dependence of the peak positions on [i is extremely weak. 



(58) 
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